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ABSTRACT 

We present results from multi-wavelength simultaneous X-ray and radio observations of the 
black hole X-ray binary V404 Cyg in quiescence. Our coverage with NuSTAR provides the 
very first opportunity to study the X-ray spectrum of V404 Cyg at energies above 10 keV. The 
unabsorbed broad-band (0.3-30 keV) quiescent luminosity of the source is 8.9 xlO 32 erg s -1 for a 
distance of 2.4 kpc. The source shows clear variability on short time scales (an hour to a couple of 
hours) in radio, soft X-ray and hard X-ray bands in the form of multiple flares. The broad-band 
X-ray spectra obtained from XMM-Newton and NuSTAR can be characterized with a power-law 
model having photon index T=2.12±0.07 (90% confidence errors); however, residuals at high 
energies indicate spectral curvature significant at a 3cr confidence level with e-folding energy of 
the cutoff to be 2C% 20 keV. Such curvature can be explained using synchrotron emission from the 
base of a jet outflow. Radio observations using the VLA reveal that the spectral index evolves on 
very fast time-scales (as short as 10 min.), switching between optically thick and thin synchrotron 
emission, possibly due to instabilities in the compact jet or stochastic instabilities in accretion 
rate. We explore different scenarios to explain this very fast variability. 

Subject headings: black hole X-ray binaries: general — low mass X-ray binary: individual (V404 Cygni) 
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1. Introduction 

Low-mass X-ray binaries (LMXBs) are systems 
containing a compact object, either a black hole 
(BH) or a neutron star (NS), accreting material 
from a low mass (less than a solar mass) late- 
type companion star. These systems often show 
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transient behavior with long intervals (years to 
decades) of low luminosity emission (namely the 
quiescent state) between bright outburst states, 
lasting from weeks to months. X-ray luminosities 
can typically have values in the range of 10 31-33 
erg s -1 during quiescence, dramatically increasing 
by several orders of magnitude during outbursts. 
A disk instability model, originally proposed to 
explain dwarf nova outb ursts, can broadly ex¬ 
plain outbursts in LMXBs ([Cannizzoll 19931 : lLasota 
2000l ). Several models have been proposed to ex- 


by Swift/BAT on 2015 June 15 ([Barthelmv et al 
20151 ). Following the Swift BAT trigger, the 


plain the origin of the quiescent X-ray emission in 
these systems. During the low-luminosity phase 
the accretion rate from the companion drops be¬ 
low a critical value, resulting in drastic structural 
changes in the accretion disk. The inner part 
of the disk can evaporate forming an advection- 
dominated accretion flow (ADAF; iNaravan fe Yi 
Il994l ) that is radiatively inefficient. Another pos¬ 
sibility invokes a jet-dominated phase where a sig¬ 
nificant amount of energy is released in an out ¬ 
flow instead of in an ADAF ([Fender et al.l 12003 ) 


lYuan et ah! ( 20051 ) suggested another interesting 
possibility of a hybrid of jet and ADAF models 
to explain spectral and timing cha racteri stics of 


XTE J1118+480. Blandford fe Begelman (jl999f ) 


proposed a modified version of ADAF in which 
powerful winds carry away mass, angular momen¬ 
tum and energy from the accreting gas. Typically, 
it has been presumed that accretion onto a qui¬ 
escent BH happens via some kind of radiatively 
inefficient accretion flow where X-rays are emit¬ 
ted via a population of hot electrons. However, it 
is still not clear if the quiescent BH always pro¬ 
duces steady collimated jets. Owing to the quies¬ 
cent nature of these BH systems, it is difficult to 
obtain good quality multi-wavelength data, with 
V404 Cyg being one of the few sources with a high- 
quality broad-band SED in the quiescent state. 

V404 Cyg has been used as a primary test 
bed for many of the above models owing to its 
close proximity and high luminosity durin g quies¬ 
cence ([Naravanll 19961: INaravan et al.lll997l ). X-ray 
emission was discovered from V404 Cyg during 
its outburst in 1989 by the Ginga X-ray satel¬ 
lite ( Makinol 19891 ). Before that, it had under¬ 
gone two outbursts, in 1938 and 1956, which were 
only seen in optical due to lack of space instru¬ 
ments. Recently, the source has been observed in 
outburst for a fourth time when it was detected 


source has been extensively monitored at multi¬ 
wavelengths in cluding optical, rad i o, UV and X- 
rays ( see e.g. Gazeas et al.l 12015 . Gandhi et al. 
120151 iKuulkers 20151 and several other related 


ATels.). V404 Cyg is believed to contain a K0(±1) 
III- V secondary star orb iting around a ~10 M 0 
BH ( Shahbaz et al.lfl996t ) with an orbital period of 
6.47 days. More recently. Ixhargharia et al . ( 201 0) 
estimated a slightly lower value of 91 q' 6 M© for 
the BH using near-infrared spectroscopy. At a dis¬ 
tance of 2.39T0. 14 k pc, obtained using parallax 
measurement ( Miller-Jones et al1l2009l ). V404 Cyg 
contains the most luminous quiescent BH, brighter 
by ab out two orders of magnitude ( Garcia et al.l 
2001 ). Typically, V404 Cyg is known to vary 


in X-rays by a factor of 2- 10 on time scales of 
a few hour s in q u iescence (I Wagner et al. 1994 
Kong et al. _ 2002 : _ Bernardini fe Cackett 20141 ). 


However, Hvnes et al. ( 2001 ) reported X-ray flux 


variation with unusually large amplitude (factor of 
>20 in 60 ks). The X-ray variations are well cor¬ 
related with that observed in the optical H a line, 
suggesting that the optical variations are powered 
by X-ray irradiation. Strong variability has also 
been reported at radi o waveleng ths, represented 
as large flares and dips ( Hvnes et al .1120091 ) with no 
clear correlations between X-ray and radio emis¬ 
sion levels. 

In this paper we present results from our recent 
X-ray and radio observations of V404 Cyg in the 
quiescent state. The next section describes the 
observations and data analysis for XMM-Newton , 
NuSTAR and the VLA. Sections 3 and 4 present 
results from timing and spectral analyses. We dis¬ 
cuss results in Section 5 and summarize them in 
Section 6. 

2. Observations and Data Reduction 

We observed V404 Cyg with XMM-Newton , 
NuSTAR (Nuclear Spectroscopic Telescope Array) 
and the VLA (Very Large Array) as part of a 
multi-wavelength campaign to study its broa d¬ 
band properties in quiescence. NuSTAR ([Harrison et al 
2013) has two coaligned hard X-ray optics mod¬ 


ules with corresponding focal planes, referred to 
as FPMA and FPMB, that cover 3-79 keV energy 
range. XMM-Newton employs X-ray optics with 
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PN and MOS CCD detectors that cover the 0.3-10 
keV energy range. The VLA provides coverage in 
the radio part of the electromagnetic spectrum in 
4.74-7.96 GHz. V404 Cyg was observed simulta¬ 
neously with NuSTAR and XMM-Newton on 2013 
October 13-14 for about 150 ks and 52 ks, re¬ 
spectively. According to our original observation 
strategy, the VLA observations were supposed 
to be simultaneous with the XMM-Newton and 
first NuSTAR observations. However, radio ob¬ 
servations were canceled due to a VLA shutdown. 
Therefore, we obtained additional simultaneous 
X-ray and radio data with NuSTAR and the VLA 
on 2013 December 2 with about 25 ks exposure 
for NuSTAR. Details of the X-ray observations are 
listed in Table [U 

2.1. XMM-Newton 

We obtained a 52.6 ks observation of V404 Cyg 
with XMM-Newton (ObsID 0723310201) on 2013 
October 13-14, covering part of the first NuSTAR 
observation (see Tabled]). We used the Scientific 
Analysis System (SAS) vl3.0.1 software package 
to reduce the data from the pn, MOS1, and MOS2 
instruments, which were all operated in full frame 
mode with the medium blocking filter. The 10-12 
keV light curve from the full pn detector shows 
that proton flares occurred at the beginning and 
end of the observation, and we filtered out the time 
when 10-12 keV pn count rate was greater than 
1 c/s, leaving a net exposure time of 31.5 ks for pn 
and 34.9 ks for each of the MOS units. The part 
of the observation used to make energy spectra 
started at 2013 October 13, 20.2 h UT and ended 
at 2013 October 14, 6.0 h UT. 

We produced event files using epproc and 
emproc, and then applied the standard event fil¬ 
tering. For pn, this means keeping events with the 
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Fig. 1.— XMM-Newton pn light curve of V404 
Cyg in the 0.5-10 keV energy band with 100 s 
time binning. 
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Fig. 2.— NuSTAR FPMA light curves for three 
epochs of observations in the 3-25 keV energy 
band with 6000 s time binning. The solid hori¬ 
zontal line in the top panel indicates the time of 
the XMM-Newton observation. 
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Fig. 3.— XMM-Newton and NuSTAR light curves 
(top panel) for the overlapping section of the ob¬ 
servation (horizontal line in Figure [2] top panel) 
with 6000 s time binning. The two count rates are 
plotted against each other (bottom panel) to see 
correlation between the two light curves. 
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Table 1 

X-RAY OBSERVATION LOG FOR V404 CYG. 


Observatory 

ObsID 

Start time 

UTC 

Stop time 

UTC 

Exposure 

(ks) 

NuSTAR (Epoch 1) 

30001010002 

2013-10-13 12:06:07 

2013-10-14 21:46:07 

61.4 

NuSTAR (Epoch 2) 

30001010003 

2013-10-14 21:46:07 

2013-10-16 23:16:07 

90.9 

NuSTAR (Epoch 3) 

30001010005 

2013-12-02 17:36:07 

2013-12-03 06:21:07 

24.9 

XMM-Newton 

0723310201 

2013-10-13 18:23:24 

2013-10-14 07:38:02 

52.6 


“FLAG” parameter set to zero, the “PATTERN” 
parameter less than or equal to 4, and applying the 
#XMMEA_EP data selection. For the MOS units, 
the filtering criterion is a “PATTERN” parameter 
less than or equal to 12 along with applying the 
#XMMEA_EM data selection. The pn and MOS 
images clearly show the presence of a point source 
at the known location of V404 Cyg, which is the 
brightest source in the field. 

We extracted a 0.5-10 keV pn light curve, a 0.5- 
lOkeV pn spectrum and 0.3-10 keV MOS spec¬ 
tra using a circular extraction region centered on 
V404 Cyg with a radius of 720 pixels (0.6 ar- 
cminutes). We used nearby source-free rectan¬ 
gular regions to estimate and subtract the con¬ 
tribution from background. We rebinned the pn 
spectrum to have at least 100 counts per bin and 
the MOS spectra to have at least 50 counts per 
bin. We produced pn and MOS response matrices 
with rmf gen and arf gen. The average count rates 
for pn, MOS1, and MOS2 are 0.105 ± 0.002 c/s, 
0.036 ± 0.001 c/s, and 0.035 db 0.001 c/s, respec¬ 
tively. 

2.2. NuSTAR 

We reduced the NuSTAR data using the 
NuSTAR Data Analysis Software (NuSTARDAS 
vl.4.1) available as part of HEASOFT v6.16 and 
the latest CALDB files. The NuSTAR observa¬ 
tions were taken during intervals of normal Solar 
activity, and we used standard filtering to remove 
periods of high background during South Atlantic 
Anomaly (SAA) passages and Earth occultation. 
We created cleaned, calibrated event files using the 
NUPIPELINE script with standard settings. The 
source spectra and light curves were extracted 
using a 50" radius circular region around the 
source. For background estimation on FPMA, 
we used nuskybgd tool that accurately generates 


the background by properly taking into account 
the non-uniformityj)f background across the focal 
plane (see lWik et al.ll20l3 for details). 

FPMB, on the other hand, showed the pres¬ 
ence of faint stray light leaking into the field-of- 
view from a nearby source during epochs 1 and 2. 
The nuskybgd tool does not account for stray light; 
thus we did not use it for FPMB. Fortunately, for 
epoch 1, the source fell on a favorable position and 
we were able to remove stray light contamination 
by carefully selecting the background from a 100 " 
radius circular region and subtracting it from the 
source. However, for epoch 2 the source was right 
on the edge of the stray light pattern, making it 
difficult to remove the contamination. Therefore, 
to safely avoid any contamination issue, we did 
not use epoch 2 FPMB data for any further anal¬ 
ysis. FPMB epoch 3 data and all FPMA data are 
clean and free from any stray light issue and we 
use them for detailed scientific analysis. 

2.3. VLA 

V404 C yg was observed w ith the upgraded Jan¬ 
sky VLA jPerlev et al.ll201 ll project 13B-016, PI: 
S. Corbel) on 2013 December 2 for a total time 
on source of ~ 9 hours. The array was in the 
B-configuration and the observations were per¬ 
formed in the C band (4.74 — 7.96 GHz) with 
two 1024 MHz base-bands centered at 5.25 and 
7.45 GHz, with eight spectral windows of 64 x 
2 MHz channels for each baseband. 

The quasar 3C286 (J1331+3030) was used for 
the bandpass calibration and to set the ampli¬ 
tude scale. We used the very nearby compact 
source J2025+3343 (0?28 away from V404 Cyg) 
to calibrate the antenna amplitude and phase 
gains. We performed a cycle of 10 minutes on 
V404 Cyg and 36 seconds on the phase-referenced 
secondary calibrator. Two additional 3-minutes 
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scans were allocated on the non-polarized com¬ 
pact source J2355+4950 to calibrate the polariza¬ 
tion leakages. Flagging, calibration and imaging 
of the data followed standard procedures and were 
carried out within the Common Astronomy Soft¬ 
ware Application (CA SA, version 4.1.0) package 


(McM ulli n et al. 2007 ). 


We set the primary calibrator 3C286 flux scales 
using the task setjy and the physic al properties 
measured by IPerlev fe Butlerl ( 2013 ). The rou¬ 
tines bandpass and gaincal were used to solve 
for the complex bandpass, the amplitude, and 
the phase gains. Stokes parameter values of the 
linearly-polarised calibrator 3C286 were assigned 
with setjy thanks to its fractional linear polar¬ 
is ation and its polarisat ion position angle listed 
in Perlev fc Butlerl ( 2013h . We derived calibration 


solutions for the leakage terms using J2355+4950. 
The flux densities of the secondary calibrators 
were bootstrapped using the task fluxscale, 
based on the absolute calibration of 3C286. Fi¬ 
nally, we applied the calibration to V404 Cyg, 
linearly interpolating the gain solutions from the 
secondary calibrator J2025+3343, assuming that 
the two sources are close enough so that the gains 
are slowly varying between the two nearby field of 
views. 


We imaged V404 Cyg in all Stokes parameters 
using the Cotton-Schwab CLEAN algorithm for 
the deconvolution. Calibrated measurements were 
split into four distinct 512 MHz bands in order to 
characterize the spectral evolution of V404 Cyg 
over the bandpass. We made 480 x 480 pixel 2 im¬ 
ages with 0715 pixel width. At 7.94 GHz the max¬ 
imum baseline is ^270kA, corresponding to an an¬ 
gular scale of ^ 0776. We used the multi-frequency 
synthesis method with the Briggs’ weighting ro¬ 
bust parameter set to 0.5 for all I images in order 
to achieve the best compromise between angular 
resolution and signal-to-noise ratio. However, the 
U and Q images were deconvolved with natural 
weighting to maximize the S/N ratio. V404 Cyg 
is consistent with a point source in all images, so 
the flux densities were measured by fitting a point 
source in the image plane, and we add a 1% sys¬ 
tematic error as usual for the VLA in this fre¬ 
quency range. 
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Fig. 4.— V404 Cyg light curves computed over 
VLA radio frequency bands as labeled in each 
panel. The binning times correspond to 1 and 0.5 
scan duration (^10 min) for the upper and bottom 
panels respectively. 
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Fig. 5.— Radio (top panel) and X-ray (bottom 
panel) light curves from simultaneous observations 
with VLA and NuSTAR during epoch 3 with 3000 
s time bin. The X-ray light curve is in 3-10 keV 
energy band and radio in 4.74-7.96 GHz. 
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2.4. Radio Flux and Polarization Mea¬ 
surement 


V404 Cyg is detected with a mean flux density 
of 300±2 /i Jy (Stokes I over the full bandwidth) in 
good agreement with previous studies conducted 
during the quiescent state (e.g., Hvnes et al .112 0091 
iMiller- Jones et ah 2008 : Gallo et al. 20051) . The 
source is unresolved with a beamsize of 0.81 x 
0.79arcsec 2 at 7.7 GHz, consis tent wi th p revious 


high spatial resolution studies (.Miller-Jones et al 

Ml). 


No significant polarised emission is detected at 
the location of V404 Cyg, as no point source can 
be fitted on the images in the Stokes C7, Q and 
V parameters. By extracting the flux from a re¬ 
duced number of channels, we avoided smearing 
out of any polarised signal over a large frequency 
bandwidth due to Faraday rotation of the linear 
polarisation vector. Therefore, we computed the 
upper limit on its polarised emission by taking into 
account the rms of the corresponding images, i.e. 
1.51, 1.48 and 1.47/iJy for Stokes parameters Q, 
U and V, respective ly. The linear polarisation, 
defined as LP = ^Q 2 + U 2 , can then be con¬ 
strained not to exceed LP < 6.33 /iJy at a 3 a 
confidence level. Regarding the fractional polar¬ 
isation FP = 100TP//, an upper limit of 2.11% 
can be inferred. Similar amplitude measurements 
have been obtained for black hole candid ate^ in 
brighter, hard states (e.g. GX 339-4: [Corbel et ah 
2000f ). however, it is close to the detection limit. 


3. Timing Analysis 

Figure [T| shows the background subtracted X- 
ray light curve of V404 Cyg in the 0.5-10 keV 
energy band from XMM-Newton pn with 100 s 
time binning. The light curve clearly shows sig¬ 
nificant variability due to flaring activity from the 
source. The count rate shows variability by a fac¬ 
tor of 4-5 during the XMM-Newton observation 
on a time scale of a few hours. Figure [2] shows 
hard X-ray light curves in the 3-25 keV energy 
band from NuSTAR FPMA module. The three 
panels correspond to the three epochs of NuSTAR 
observations, with a solid horizontal line in the top 
panel indicating the interval corresponding to the 
XMM-Newton observation. NuSTAR light curves 
are binned to 6000 s, approximately corresponding 
to an orbit of NuSTAR. The source clearly shows 



Fig. 6.— 3-10 keV X-ray flux from NuSTAR 
epoch 3 observation and radio flux (4.74-7.96 
GHz). 



Fig. 7.— NuSTAR FPMA spectra for three epochs 
of observations. Black, red and green points cor¬ 
respond to epochs 1,2, and 3, respectively. 
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significant variability in the hard X-ray band, with 
similar time scale as XMM-Newton. This suggests 
that the flaring behavior continues in the hard X- 
ray band as well. During the epoch 3 snap-shot 
observation with NuSTAR , the source was caught 
mainly in a flaring state, suggesting that it was 
showing persistent flaring behavior during the en¬ 
tire span of these observations. The three epochs 
of NuSTAR observations caught about 6-7 small 
flares from this source. Figure [3] shows X-ray light 
curves from the overlapping sections of the XMM- 
Newton and NuSTAR epoch 1 observations (top 
panel) to facilitate a direct comparison between 
the soft X-ray (0.5-10 keV) and hard X-ray (3-25 
keV) variations. The two light curves show very 
similar profile suggesting a strong correlation be¬ 
tween the two X-ray bands. We plot the XMM- 
Newton count rate and the NuSTAR count rate in 
the bottom panel of Figure [3] to look for any corre¬ 
lation between the two light curves in top panel. It 
is clearly visible that the X-rays in two bands are 
correlated. In order to quantify this correlation, 
we measured the Pearson’s correlation coefficient 
between the two light curves to be 0.9^0 *.2 (90% 
confidence errors). All the errors in this paper 
are quoted with 90% confidence limit for a sin¬ 
gle parameter of interest, unless otherwise noted. 
A value of 1 for this coefficient indicates a strong 
positive correlation and a value of 0 indicates no 
correlation. Hence, a value of 0.9 suggests a strong 
correlation between 0.5-10 keV and 3-25 keV X- 
rays, as expected. 

In order to search for shorter time scale vari¬ 
ability in radio, the calibrated VLA data were split 
into time segments with a duration chosen to be as 
small as possible for the detection to be considered 
significant. The source was fit in the image plane 
and the resulting radio light curves of V404 Cyg 
are shown in Figure IU highlighting the flux density 
variations across the observation. V404 Cyg shows 
evidence of fast radio variability with a flux den¬ 
sity ranging from 0.09 to 0.64 mJy (integrated over 
the whole frequency band). There is the sugges¬ 
tion of a small flaring activity around 100 min after 
the beginning of the observation, with a 0.64 mJy 
peak level. This flare is characterized by a fast rise 
of less than 5 min, followed by a slower decrease 
over ^30 min. 

In order to directly compare the source variabil¬ 
ity in the radio and X-ray bands, we plot their re¬ 


spective light curves in Figure [5] The light curves 
are from the simultaneous radio (VLA in 4.74-7.96 
GHz) and X-ray observations ( NuSTAR epoch 3 
in the 3-10 keV energy band) and binned to 3000 
s. The source exhibits significant variability at 
both radio and X-ray energies. The radio flux den¬ 
sity changes by a factor of ^3 during the time of 
observation. Although clear variability is present 
in both bands, there is no strong correlation be¬ 
tween radio and X-ray fluxes, partly due to the 
poor signal-to-noise ratio of the X-ray data. Once 
again we measured the Pearson’s correlation coeffi¬ 
cient between radio and X-ray time series, to quan¬ 
tify any correlation between the two. The corre¬ 
lation coefficient between radio and X-ray bands 
is 0.4lo J- A value of 0.4 suggests a weak positive 
correlation, though the large errors make it sta¬ 
tistically insignificant. Figure [6] shows X-ray and 
radio flux variation for better visualization. Thus 
we conclude that there is no compelling evidence 
for correlation between the simultaneous radio and 
X-ray data analyzed here. 


4. Spectral Analysis 


4.1. X-ray Spectral Analysis 


We studied the spectra for each X-ray observa¬ 
tion separately in order to investigate any spec¬ 
tral variation from epoch to epoch. Throughout 
this work we performed spectr al modeling using 
XSPEC vl2.8.2 (lArnaudl Il996h . Figure 0 shows 
the NuSTAR FPMA spectra for the three epochs 
of observations. It is readily visible that the 
spectra follow each other very closely, suggesting 
there is very little spectral variation among these 
three epochs. In order to quantify any spectral 
variation, we examined each spectra separately 
(three NuSTAR spectra and one XMM-Newton ). 
We started with a simple model consisting of a 
power law modified by photoelectric absorption 
due to interstellar medium. To account for the 
absorption due to intervening material, we used 
the tbab s spectral model wit h updated solar abun¬ 
dances ( Wilms et al.l 120001 ) and photoionizatio n 
cross-sections as described in lVerner et al.l ( 19961 ). 


For XMM-Newton , the EPIC (pn and MOS) 
spectra were fitted simultaneously. We obtained 
a best fit value of photon index of r=2.04±0.08 
and absorption 7Vh=(1.18±0.09) x 10 22 cm -2 . 
The absorbed 0.3-10 keV flux is 5.5 x 10 -13 


7 














Table 2 

Best fit spectral parameters from joint fit to XMM-Newton and NuSTAR data. 


Model 

N n 

r/r p 

E cu t/kT e 

Norm 

X 2 /dof 

Flux* 


(10 22 cm -2 ) 


(keV) 

l 

o 



Power-law 

1.25±0.08 

2.12±0.07 

— 

2.20±0.20 

173/200 

3.31 

Cut-off power-law 

1.14±0.10 

1.86=1=0.15 

20+ 20 

1 Q2+ 0 - 22 
l.UZ —o 20 

161/199 

3.39 

Bremsstrahlung 

0.84±0.05 

— 

6 9+ 0 ' 7 
°- y -0.6 

1 e;q+ 0 - 08 

l.OO — o Q7 

199/200 

3.60 

CompTT 

0.78±0.06 

oo 

1 + 

° ? 

oo -i 

4 0" 1 " 1 ’ 2 
^• U -0.6 

0-6(rt“;“ 

160/199 

3.42 


* Absorbed flux in 3—10 keV energy band with units of 10 13 erg cm 2 s' 



Energy (keV) 

Fig. 8.— Top panel: Broad-band V404 Cyg spec¬ 
tra from XMM-Newton (EPIC-pn in black, MOS1 
in red and MOS2 in green) and NuSTAR (FPMA 
in blue and FPMB in light blue). The solid lines 
represent best fit power-law models. The lower 
panels show the ratio of the data to different spec¬ 
tral models, as labeled. 
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which agrees well with the pre¬ 
viously reporte^_yalue_o^5.8x 10 -13 erg s -1 
cm -2 (Bernardini & Cackett 201 4). The absorp¬ 
tion obtained from XMM-Newton is a little higher 
th an the Galactic value of 8.10 x 10 21 cm -2 
( Dickev fe Lockman 1990h . While performing fits 
with individual NuSTAR spectra, we fixed TVh to 
the value obtained with XMM-Newton as NuS- 
TARA energy range is not sensitive to such ab¬ 
sorption values. The best fit values for the power- 
law index for the three NuSTAR observations are 
2.542.28±0.17 and 1.90±0.24, respectively. 
The spectrum apparently is battening going from 
epoch 1 to epoch 3, however, the individual er¬ 
rors on the power-law indices overlap, suggesting 
that the difference may not be statistically signifi¬ 
cant. In order to investigate this we fitted the data 
from each epoch simultaneously with an absorbed 
power-law model. Ah the parameters were linked 
between epochs and only the overall normalization 
constants were allowed to vary. This also provided 
an excellent fit to the data, with a reduced y 2 of 
1.07 for 95 degrees of freedom and a best fit power- 
law index of r=2.35±0.12. This suggests that any 
variations in the spectral parameters during the 
three epochs are not significant and the spectra 
can formally be considered to be consistent with 
each other. Further, we have also derived hardness 
ratios between the two NuSTAR energy bands, 3- 
8 keV (soft) and 8-30 keV (hard) to test for any 
X-ray spectral variations. We found no variation 
in the hardness ratio (hard/soft), indicating that 
there are no significant spectral changes in the X- 
ray data. This is consisten t with the previou s X- 
ray stud ies u si ng Chandra ( Hynes ct al.l2009j ) and 
Swift 


Bernardini fe Cackettll2014l ). 

In order to evaluate the broad-band X-ray char¬ 
acteristics of V404 Cyg during quiescence, we si- 































multaneously fitted XMM-Newton and NuSTAR 
data to achieve energy coverage from 0.3-30 keV. 
As shown above, there is no significant spectral 
variation between the three epochs of NuSTAR 
observations, therefore we combined the data from 
the three NuSTAR epochs to improve the signal- 
to-noise ratio. We used the ftool addascaspec to 
combine the NuSTAR data, which properly com¬ 
bines the corresponding background and response 
files. In order to account for the cross-calibration 
between different instruments, we included a con¬ 
stant multiplicative factor in all the spectral mod¬ 
eling described below. This constant factor was 
held fixed to 1.0 for XMM-Newton EPIC-pn and 
allowed to vary freely for other instruments. The 
EPIC-pn to NuSTAR normalization diff erenc es 


are b elow the expected <10% level (Madsen et al 
2fii rfi . 


We employed several different spectral mod¬ 
els to characterize the broad-band spectra of 
V404 Cyg during quiescence, including an ab¬ 
sorbed power-law model, a cutoff power-law 
model, a thermal bremsstrahlung model and a 
thermally comptonized continuum model (see 
Figure [8|). The best fit spectral parameters ob¬ 
tained through joint fitting of XMM-Newton and 
NuSTAR data are listed in Table [2] All these 
models provide statistically acceptable fits to the 
data. The power-law model provides a best fit 
value of the photon index of T=2.12±0.07. The 
unabsorbed flux in 0.3-30 keV energy band is 
(1.29±0.07) x 10 -12 erg cm -2 s -1 and the corre¬ 
sponding luminosity is (8.9±0 .5) x 10^_ erg s" 1 
for a distance of 2.4 kpc ( Miller-Jones et al1l2009h . 
Although the power-law model provides an accept¬ 
able fit (x 2 =0.86), the residuals show a hint of a 
roll-over in the high energy data above ^10 keV. 
We therefore tried a power-law model with a high- 
energy exponential cut-off (cutoffpi in XSPEC) 
to test if it accounted for the high-energy residu¬ 
als. This model, with an e-folding energy of the 
cut-off of 20^7° keV, further improves the fit, with 
Ay 2 =12 for one additional degree of freedom. The 
cut-off energy, being very close to the high end of 
the spectral data, is not tightly constrained. 

In order to assess the significance of the poten¬ 
tial cutoff, we performed a series of spectral sim¬ 
ulations. Using the same responses, background 
files and adopting the same exposure times as 
the data used here, we simulated 10,000 sets of 


XMM-Newton (pn, MOS1, MOS2) and NuSTAR 
(FPMA, FPMB) spectra with the FAKEIT com¬ 
mand in XSPEC. These simulations were based 
on the best-fit power-law continuum, and each of 
the simulated datasets was rebinned in the same 
manner and analyzed over the same bandpass as 
adopted for the real data. We then fit each of 
the combined datasets with an absorbed power- 
law continuum with and without an exponential 
cutoff, again in the same manner as the real data, 
and noted the improvement in fit provided by the 
former over the latter. Of the 10,000 datasets 
simulated, only 28 showed a chance improvement 
equivalent to or greater than that observed, imply¬ 
ing that the detection significance of the cutoff in 
the real data is at the 3cr level. Next, we performed 
another set of simulations with the best fit cutoff 
power-law model as input, in order to assess the 
expected distribution of Ay 2 when compared to a 
simple power-law model given the available data. 
We found that it peaks near the observed value 
and this holds true even if we take into account 
that the power-law model already provides a good 
fit with reduced y 2 <0.9. All of these simulations 
favor the detection of a cutoff with 3 a significance 
in our broad-band XMM-NewtonTNuSTAR data. 
Such a feature has been measured for the first time 
in V404 Cyg during quiescence. 

We next used a physically motivated thermal 
bremsstrahlung model to characterize the broad¬ 
band spectra. Again, this model provides a sta¬ 
tistically acceptable fit to the data. The derived 
plasma temperature is 6.9 <q 6 keV. However, this 
fit is poorer by Ay 2 =33 compared to the cut¬ 
off power-law model. Thus the cut-off power- 
law is the favored spectral model for these data 
among the above three models. In addition, we 
also tried a thermally Comptonized continuum 


Table 3: Summary of V404 Cyg Stokes I measure¬ 
ments. A point source is fitted to each frequency- 
dependent image integrated over a 512 MHz band¬ 
width,_ 


Frequency 

(GHz) 

Flux density 

(£iJy) 

Beam size 
(arcsec 2 ) 

5 

314.7 ±6.2 

1.25 x 1.21 

5.5 

307.2 ±5.5 

1.13 x 1.09 

7.2 

282.8 ±4.0 

0.84 x 0.83 

7.7 

281.8 ±4.0 

0.81 x 0.79 
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Fig. 9.— V404 Cyg radio spectrum as described 
in Table 02 The dotted black line and the gray 
shaded area represent the fitted power-law shape 
spectrum and the 95 % confidence interval respec¬ 
tively. The spectral index for the power-law is 
<a=-0.27=b0.03. 



Fig. 10.— Variation of the spectral index a over 
the length of the observation obtained by fitting a 
power-law (f u oc i/*) to the data split in ~ 10 min 
time interval. 


model, comptt , that describes Comptonization of 
soft photons in a hot plasma from a Comptoniz- 
ing corona a nd tiny i nner region of an ADAF flow 
dTitarchuk fe Lvnbarskiil Il995l : iHua fc Titarchukl 


19951 ). The comptt model has a higher number 


of free parameters, and not all of them are con¬ 
strained simultaneously. The seed photon temper¬ 
ature kTo varies between 0 to 0.4 keV with a best 
fit value of 0.32 keV, which is in excellent agree- 
me nt wit h the previous measurement of 0.33 keV 
bv iBernardini fe Cackett ( 2014 ). Hence we fixed 
the seed photon temperature to 0.32 keV during 
the fit. As listed in Table [2j the comptt model is 
characterized by a coronal electron temperature of 
~4 keV and an optical depth of ~5. The goodness 
of fit with this Comptonization model is very sim¬ 
ilar to that found for the simple cutoff power-law 
model (see Table [2j) and shows very similar resid¬ 
uals. Therefore, we have not shown it in Figure [8] 
for the sake of clarity. 

Visual inspection of the X-ray spectra indi¬ 
cates that the data are featureless, showing only 
a smooth continuum over the broad X-ray band. 
However, there are predictions for the presence 
of X-ray emission li nes from the h ot plasma i n 
an ADAF-like flow dNaravan fe Raymond! Ill 9991) . 
In order to derive an upper limit on the highly 
ionized Fe XXV K a line at 6.7 keV, we added a 
Gaussian component to our best fit cutoff power- 
law model. We obtained an upper limit on the 
equivalent width of ^200 eV (0.1 keV fixed width 
and 90% confidence limit) from the broad-band 
XMM-Newton and NuSTAR spectra. This is 
so mewhat high er than the upper limit derived 
bv lBradlev et al.l ( 20071) . However, it is nearly as 
large as the prediction of ^230 eV for the sum 
of the equivalent widths for t h e tota l grou p of Fe 
emissoin lines ( Naravan fe Raymond! 11999 ). 


4.2. Radio Spectral Analysis 

In order to characterize the spectral energy dis¬ 
tribution (SED) at radio energies, we measured 
the flux density of V404 Cyg through four distinct 
frequency bands over the whole observing period. 
We report the results in Table [3j We use this in¬ 
formation to constrain the radio spectral shape by 
fitting the function f v oc z/*, where f v is the flux 
density and v is the frequency. Figure [9] shows 
the radio spectrum (black points) fitted with a 
power-law model (dotted line) and the correspond- 
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ing 95% confidence region (gray shaded area). We 
obtain a spectral index of a = —0.27 =b 0.03 (lcr 
error) for the average spectrum across the whole 
observation. The same fitting procedure was ap¬ 
plied to two frequency bands to derive the evolu¬ 
tion of the spectral shape during the observation. 
Figure HOl shows the spectral index variation over 
^10 min time bins. Interestingly, the spectral in¬ 
dex changes rapidly between optically thick (0.0 to 
0.5) and thin (-1.0 to -0.5) synchrotron emission 
regimes. In order to make sure that such extreme 
variability is associated with V404 Cyg, we per¬ 
formed a similar analysis on nearby field sources. 
We found that the flux density and the derived 
radio spectral index of the field sources were con¬ 
stant over the course of the observation, confirm¬ 
ing that the observed spectral variability displayed 
in Figure [10] is associated with V404 Cyg. 

5. Discussion 

We obtained simultaneous multi-wavelength 
data for V404 Cyg during its quiescent state, al¬ 
lowing us to extend such studies beyond 10 keV 
for the first time. 


5.1. Comparison with Previous Studies 


Over the course of our observations, V404 Cyg 
showed strong variability at soft X-ray, hard X- 
ray and radio wavelengths, represented as multiple 
flares with each flare lasting for a few hours. Sim¬ 
ilar temporal variability has been observed previ¬ 
ously in long term mo nito ring of_this source with 
Swift XRT ( Bernardini &; Cackettl 1201 4) in the 
0.3-10 keV energy band. That study showed that 
V404 Cyg is highly variable on a wide range of 
time-scales from tens of minutes to years in 34 
Swift obser vations of the source in quiescence. 
Hvnes et al.l ( 2004 ) also reported dramatic X-ray 
flux variability with much higher amplitude (more 
than a factor of 20) using 2003 Chandra obser¬ 
vations in the 0.3-7 keV band. Our sensitive 
NuSTAR observations clearly demonstrated that 
the flaring behavior of V404 Cyg continues be¬ 
yond 10 keV. In addition, new VLA radio light 
curve, taken simultaneously with the NuSTAR 
hard X-ray data, also showed strong variability 
on short tim e scal es, similar to that reported by 
Hynes et al. ( 2009 ). However, unlike these previ¬ 
ous radio light curves, our new radio data show a 


persistent emission component with the flux level 
never going to zero. Although the simultaneous 
radio and X-ray data show strong variability, we 
do not find any significant correlation (or anti¬ 
correlation) between the two data sets. This is 
again c onsistent with pre vious multi-wavelength 


2009h . 


studies (jHvnes et al 

The quiescent X-ray spectra of V404 Cyg 
have been studied previously with XMM-Newton ., 
Chandra and Suz aku (see Bradley et al. _ 20071 


Corbel et al. 20081 and Reynolds et al.l 12014 ) and 


are well characterized using absorbed power- 
law models. We measured the power-law index 
r=2.12±0.07 using our broad-band XMM-Newton 
and NuSTAR data. This measurement is in excel¬ 
lent agreement with previously reported power- 
law spectra from Chandra (r=2.17±0.13) and 
XMM -Newton (r=2.09±0.08). iRnvnolds et al. 
( 2014 ) recently studied quiescent X-ray spec¬ 
tra for a sample of accreting BHs, including 
V404 Cyg. They reported best fit power-law 
index of 2.05T0.06 from a simultaneous fit to 
XMM-Newton , Chandra and Suzaku observations, 
which agrees with our measurement within 90% 
errors. The broad-band XMM-Newton and NuS¬ 
TAR spectra, along with excellent sensitivity, 
show a presence of cutoff at around ~20 keV with 
more than ~3cr significance. A similar cutoff (at 
around 10 keV) was recently detected by NuSTAR 
in the quiescent spectr a of the low ma ss X -ray bi¬ 
nary Cen X-4 ( Chakrabartv et all 2014 ) which 
these authors attribute to hard X-ray emission 
arising either from a hot, optically thin corona or 
from hot electrons in an optically thin radiatively 
inefficient accretion flow (RIAF). The broad-band 
spectrum of Cen X-4 is well characterized by a 
thermal bremsstrahlung model with kT e = 18 keV. 
For V404 Cyg, the bremsstrahlung model also 
provides a statistically acceptable fit, though, the 
best fit temperature, kT e = 6.91q 6 keV, is sig¬ 
nificantly l ower t han that obse rve d in Cen X-4. 
Recently, llBernardini fc Cackettl ( 2014 ) analyzed 
Swift XRT spectra of V404 Cyg and reported that 
fits with a bremsstrahlung model find the temper¬ 
ature between 3.3 to 4.7 keV as a function of count 
rate, which is sl ightly lower than our measure¬ 
ment. However, iRevnolds et al. ( 2014 ) reported 
a bremsstrahlung temperature of 6.16lggg keV 
using Suzaku XIS data, in agreement with our 
measurement. 
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Corbel et al.l (112008) revisited the radio-X-ray 
flux correlation for V404 Cyg by re-analyzing 
some of the 1989 outburst data and adding two 
points from their 2003 Chandra/Y LA data and 
2000 Ginga/Y LA observations during quiescence. 
They found a strong correlation between the two 
energy bands and conclude that the correlation 
observed during the decay of the 1989 outburst 
holds down to the quiescent state. Using our 
new simultaneous NuSTAR/YLA observations, we 
were able to add one more point to that corre¬ 
lation, shown jin Figure [TlJ Earlier data were 
taken from Corbel et al. (12008[ ) (see their Figure 4) 
and we did not attempt any fitting. Our new 
data point lies right on their best fit curve, nicely 
following the long-term correlation from previ¬ 
ous data. The quiescent flux level of V404 Cyg 
during our observations is higher compared to 
their Chandra/Y LA observation, but it is signifi¬ 
cantly lower than the Ginga/Y LA measurements. 
This long-term correlation between X-ray and ra¬ 
dio flux is indicative of a common or related ori¬ 
gin for the two components. Although the radio 
and X-ray fluxes are correlated on long timescale, 
the short timescale variations during the observa¬ 
tions become more complicated and any correla¬ 
tion mny_^e_smoothed_out by internal jet variabil¬ 
ity ( Gleissner et al.l 20041) . 

The quiescent X-ray flux of V404 Cyg shows sig¬ 
nificant variability on a time scale of a few years, 
as reported in previous studies. iRevnolds et al 


(120141 ) reported that the 0.3-10 keV unabsorbed 


X-ray flux varies in the range of (0.8-3.42) x 
10“ 12 erg s -1 cm -2 using various observations 
between 2000-2009 from Chandra , XMM-Newton , 
and Suzaku. We measured the unabsorbed 0.3-10 
keV flux to be 1.1 x 10“ 12 erg s -1 cm -2 which is 
well within the above range, however it is slightly 
lower compared to the average flux of 1 .7 x_10 -12 
erg s -1 cm -2 , reported bv iBernardini fe Cackettl 
(HI using 75 d ays of Swift monitor ing obser¬ 


vations in 2012. Naravan et al. (Il997h reported 


unabsorbed 1-10 keV flux of 8.2 x 10 -13 erg s“ 
cm -2 using 1994 ASCA observations. Thus, the 
X-ray flux of V404 Cyg shows significant variabil¬ 
ity; however, its spectral shape remains fairly sta¬ 
ble with T=2.1-2.3. 


5.2. Hard X-ray Emission Mechanisms 

We can use the measured spectral cutoff to 
make inferences about the nature of the accretion 
flow onto the bl ack h ole. We follow thej^eneral ap¬ 


proach used by Chakrabartv et al. (2014) in inter¬ 


preting the spectral cutoff in the quiescent accret¬ 
ing neutron star Cen X-4. A key difference here is 
that the observed X-ray luminosity in black hole 
systems provides only a lower limit on the mass 
accretion rate, since matter can be accreted onto 
the black hole without radiating. In the discussion 
below, we assume a black hole mass of 10 M© and 
a distance of 2.4 kpc, and we write the Eddington 
luminosity and accretion rate for a 10 M© black 
hole as Le = O.lM^c 2 = 1 x 10 39 erg s -1 . The 
observed X-ray luminosity is thus L x ~ 10 Le- 
The hard X-ray emission in BH binaries is of¬ 
ten ascribed to a Comptonizing corona of hot elec¬ 
trons. Our comptt fit yielded an electron temper¬ 
ature of kT e = 4.0 keV and an electron scattering 
optical depth of r = 4.8. This is a fairly high opti¬ 
cal depth, and it is related to the electron density 
n e by 

T = (? T J n e (r)dr , (1) 

where gt is the Thomson cross-section. We as¬ 
sume a quasi-spherical, RIAF inside a transition 
radius at r t = lO 4 ^, where Rs = 2 GM/c 2 = 
30 km for a 10 M© black hole. We then find that 
an accretion rate of M « 0. 02Me is required to 
produce the measur ed value of r. (We used Equa¬ 
tions (5) and (8) of IChakrabartv et al.l ( 2014 ) to 
relate n e and M, taking g = 0.1 and g = 0.6.) 
The resulting electron population has a high scat¬ 
tering optical depth but is optically thin to free- 
free absorption, and will thus itself be a source 
of bremsstrahlung emission. At a temperature of 
kT = 3.9 keV, this would produce an X-ray lumi¬ 
nosity 60 times brighter than what we observe. If 
we allow for the possibility of an outflow by setting 


M(r) = ( - ) M t 


( 2 ) 


with p > 0 ( Blandford fc Begelmanll19991 ), the dis¬ 
crepancy is even larger. We conclude that ther¬ 
mal Comptonization cannot be responsible for the 
spectral cutoff. 

Thermal bremsstrahlung emission from the ac¬ 
cretion flow is, however, a viable possibility. Our 
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bremsstrahlung fit has an emission measure 


/ 


l dV = 4.7 x 10 57 cm” 3 


(3) 


for a distance of 2.4 kpc. For the p = 0 case 
in Equation (2), an accretion rate of M ~ 
10 ~ 3 M E is s ufficient to yield this emission mea¬ 
sure. iMenou et al.l ( 119991 ) suggested that V404 


Cyg during its quiescent state is basically in the 
RIAF regime when mass loss via wind is included. 
If we allow for an outflow (p > 0), then we require 
M t ~ 10 ~ 2 Me at the transition radius. These 
are plausible accret ion rates for the RIAF regime 
( Menon et al.l[l999l ). 


Another possibility is synchrotron and/or syn¬ 
chrotron self- Compton emission fr om the base of 
a jet outflow (Ma rkoff et al.ll2~005l) . Such models 
are competitive with Compton coronal models for 
the X-ray emission in black hole binaries and are 
capable of producing curvature in the power-law 
spectrum which could be consistent with our ob¬ 
served cutoff. Detailed modeling will be necessary 
to evaluate the viability of this model in V404 Cyg. 

5.3. Rapid Radio Spectral Variation 

The hundreds of seconds to ks time-scale flar¬ 
ing behavior of V404 Cyg in quiescence has al¬ 
ready been reported at several wavele ngths, in¬ 
cluding at 8.4 GHz in the radio dMiller-Jones et al 
2008 : iHvnes et al, 


2009h . Our new sensitive ra¬ 


dio observations with the VLA confirm this be¬ 
havior, but they also allow us for the first time to 
have additional radio spectral information during 
these periods of weak activity. The spectral in¬ 
dex (Figure fT0|) evolved on a very fast time-scale 
(as short as 10 minutes) with clear switching be¬ 
tween regimes of optically thick and optically thin 
synchrotron emission. 

In order to explain the observed spectral vari¬ 
ability, we explored several possibilities, one of 
them being intervening material along the line 
of sight. V404 Cyg lies behind the Cygnus 
super-bubble, composed of interstellar structures 
shaped by supernova explosions and stellar winds 
of massive stars, all estimated to be located 
bet ween 0.5 and 2.5 kpc from the Eart h 


e.g. Uvamker et al. 2001 : lAckermann et al~ 2011 ) 


Fluctuation of the electron density in these clouds 
can modify the refractive index of the plasma, in¬ 
ducing scintillation of the observed radio emission. 


However, Miller-Jones et al. ( 2008 ) estimated that 
the time-scale of refractive scintillation should be 
between 140 and 1000 hrs, given the source size 
constraints they derived from VLBI observations. 
The flux variability seen in our VLA observations 
on few minutes time-scale must therefore be in¬ 
trinsic to V404 Cyg. 

In the hard state, self-absorbed compact jets 
are usually observed with a radio spectrum that 


is ei ther flat or slightly inverted (e.g. ICorbel et al 


2000t ). However, it has been found that, when the 
compact jets are building up during the soft to 
hard state transition, they are initially con sistent 
with optically thin emiss ion (e.g., ICorbel et al. 


2013 : Kalemci et al.l 120131 ). The quiescent state 
in V404 Cyg may be similar to such a situation 
with the accretion flow not being able to sustain 
stable compact jets due either to insufficient parti¬ 
cle density and/or inefficient particle acceleration. 
One could possibly interpret the rapid spectral in¬ 
dex changes in Figure [10] as jet ignition instabili¬ 
ties, rapidly moving between development/fading 
(a < 0) and activated states (corresponding to 
flat/inverted spectra, a > 0). A dedicated cam¬ 
paign with simultaneous observations in infrared 
and radio should provide an important diagnostic 
for such behavior. 

Alternatively, the observed fast spectral modu¬ 
lations of the radio spectrum may originate from 
stochastic instabilities in the accretion flow. The 
synchrotron turnover frequency, z/*,, of a compact 
jet (usuall y at infrared frequencie s in the hard 

Russell et al 


jet (usuall y at mlra red ireq uencie 
state , e.g. ICorbel fe Fenderl 2002 


20061 : IHvnes et~ah 12009 : ICorbel et al.l 2013^ could 
move to lower fre quencies as the mass accretion 
rate, m, decreases (|Falcke et al ■ 201)4: IChatv et al 
20 lll l. In quiescence, the accretion rate ir 


V404 Cyg may fluctuate on very short time-scales. 
Rapid modulations of the magnetic field inten¬ 
sity and/or the size of the particle acceleration 
zone may cause v\, to shift over to the range 


of our observ able bandwidth (Ch atv et al.l 12011 


Gandhi et al thus producing the switches 

between optically thick and thin synchrotron emis¬ 
sion, as observed at infrared wavelen gths o n sim¬ 
ilar time-scales for other sources (IRahoui et al 


20121 : iGandhi et al.l 1201 If ). The fact that the new 
observations lie on the standard radio and X-ray 
flux correlation observed in the hard and qui¬ 
escence states probably indicates an interpreta- 
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ti on related to t he co mpac t jets. As observe d 
in iGandhi et al.l ( 2011 ) and iRahoui et al.l ( 2012h . 
even if the compact jets are strongly variable on 
short time-scales, those variations smooth out on 
longer time-scales. 

The observed radio emission may also be un¬ 
related to compact jets, which might be absent 
in quiescence. The radio flaring behavior might 
instead be related to synchrotron bubble ejection 
events, possibly triggered by magnetic reconnec¬ 
tion events above the accretion disk, or be related 
to variations in the accretion rate. In the syn¬ 
chrotron bubble scenario, the high-frequency peak 
should lead the lower low-frequency flux maximum 
( van der Laanlll96(l iHiellming fc Johnstonlll988l) . 


Considering averaging effects, this is consistent 
with the rapid switch between the observed posi¬ 
tive and negative spectral indices (see Figure H0|) . 
However, in such a case, one would expect the 
higher frequency light-curve to have higher max¬ 
ima than observed at lower frequency. This is in 
contrast to what we observe (Figure 0]), and indi¬ 
cates that the synchrotron bubble scenario is not 
the favored interpretation. 


6. Summary 


V404 Cyg shows strong variability on time- 
scales of hours to a couple of hours, at both 
X-ray and radio energies. Several small flares 
were detected at X-ray energies by both XMM- 
Newton and NuSTAR , during which the count rate 
changed by a factor of 4-5 in a few hours. On 
the other hand, the radio flares are comparatively 
short lived and asymmetric, with fast rises chang¬ 
ing flux density by a factor of ~ 7 in 5 min followed 
by a slower decay. 


Our broad-band X-ray spectral analysis using 
absorbed power-law models provides consistent re¬ 
sults with previous studies using just X-ray data 
below 10 keV. However, with NuSTAR extending 
the spectral coverage above 10 keV, we were able 
to measure a cut-off at 20+7° keV, which can be 
explained using thermal bremsstrahlung emission 
from the accretion flow or synchrotron emission 
from the base of a jet outflow. V404 Cyg is the 
second quiescent LMXB showing such a feature, 
in ad ditio n to the neutron star system Cen X-4 
( Chakrabartv et al . 2014 ). 

Our NuSTAR and VLA observations are in 


very good agreement with the long-term correla¬ 
tion between X-ray and radi o flux measured ove r 


many orders of magnitude (Corbel et ah 2008 ) 


This indicates that the origins of these two com¬ 
ponents are not completely independent on long 
timescales. 


Our new sensitive radio observations discovered 
extreme variability in the spectral index on time- 
scales as short as 10 minutes, with clear switching 
between optically thick and thin synchrotron emis¬ 
sion. Such variability could possibly be attributed 
to instabilities in the compact jet or stochastic 
instabilities in the accretion rate where the syn¬ 
chrotron turnover frequency could shift to lower 
energies during quiescence. 
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Fig. 11.— X-ray flux versus Radio flux d ensi ty 
variation reproduced from [Corbel et al. (I2C 


with addition of a new point from our recent NuS- 
TAR and VLA observations shown as a square 
box. Solid circles represent data from 1989 out¬ 
burst of V404 Cyg. 


This 2-column preprint was prepared with the AAS IAT^X 
macros v5.2. 
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